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The quantum emulation of spin-momentum coupling (SMC), a crucial ingredient for the emer¬ 
gence of topological phases, is currently drawing considerable interest. In previous quantum gas 
experiments, typically two atomic hyperfine states were chosen as pseudospins. Here, we report 
the observation of a new kind of SMC achieved by loading a Bose-Einstein condensate (BEC) into 
periodically driven optical lattices. The s- and p-bands of a static lattice, which act as pseudospins, 
are coupled through an additional moving lattice which induces a momentum dependent coupling 
between the two pseudospins, resulting in s-p hybrid Floquet-Bloch bands. We investigate the band 
structures by measuring the quasimomentum of the BEC for different velocities and strengths of the 
moving lattice and compare our measurements to theoretical predictions. The realization of SMC 
with lattice bands as pseudospins paves the way for engineering novel quantum matter using hybrid 
orbital bands. 


Spin-momentum coupling (SMC), commonly called 
spin-orbit coupling, is a crucial ingredient for many im¬ 
portant condensed matter phenomena such as topologi¬ 
cal insulator physics, topological superconductivity, spin 
Hall effects, etc [1-3]. In this context, the recent ex¬ 
perimental realization of SMC in ultracold atomic gases 
provides a powerful platform for engineering many inter¬ 
esting and novel quantum phases [4-9]. In typical ex¬ 
periments, two atomic hyperfine states act as two pseu¬ 
dospins which are coupled to the momentum of the atoms 
through stimulated Raman transitions [10, 11]. However, 
ultracold atoms in optical lattice potentials possess other 
types of degrees of freedom which can also be used to de¬ 
fine pseudospins [12, 13]. A natural and important ques¬ 
tion is whether such new types of pseudospins can be 
employed to generate SMC. 

In optical lattices filled with ultracold atoms, s- and 
p-orbital bands are separated by a large energy gap and 
can be defined as two pseudospin states. One signifi¬ 
cant difference between hyperfine state pseudospins and 
lattice band pseudospins lies in the energy dispersion of 
“spin-up” and “spin-down” orientations: the dispersion 
relations are the same for hyperfine state pseudospins, 
while they are inverted for lattice band pseudospins. It 
is well known from topological insulators and supercon¬ 
ductor physics that inverted band dispersions, together 
with SMC, play a central role for topological properties of 
materials [14-16]. Therefore, it is natural to expect that 
the inverted band pseudospins, when coupled with the 
lattice momentum, may lead to interesting topological 
phenomena in cold atomic optical lattices. Recent exper¬ 
iments with shaken optical lattices (i.e. lattices in which 
the lattice sites are periodically shifted back and forth in 
time [17]) have realized a simple coupling (Ula x coupling, 
where H is the coupling strength and a x a Pauli matrix) 
between s- and p-band pseudospins, analogous to Rabi 
coupling between two regular spins [18]. However, for the 
exploration of exotic phenomena in optical lattice sys¬ 
tems, such as Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 


phases [19, 20] and Majorana fermions [16], SMC with s- 
and p-bands pseudospins is highly desirable [21-24]. 

In our experiments we realize such s-p band SMC for 
a Bose-Einstein condensate (BEC) using a weak moving 
lattice to generate Raman coupling between s- and p- 
band pseudospins of a static lattice [25]. The moving 
lattice acts as a periodic driving field [26-31] and has 
previously been used to generate an effective magnetic 
field in the lowest 5-band of a tilted optical lattice [32, 33] . 
In our experiment, the driving frequency of the moving 
lattice is chosen close to the energy gap between s- and 
p-bands at zero quasimomentum, leading to a series of 
hybrid s-p Floquet-Bloch (FB) band structures. FB band 
structures in optical lattices give rise to interesting and 
important phenomena in cold atoms and solids [34, 35], 
as is evidenced by the recent experimental realization 
of a topological Haldane model in a shaken honeycomb 
optical lattice [36] and the observation of FB states on 
the surface of a topological insulator [37]. 

Here we show that the moving lattice generates two 
types of coupling between s- and p-band pseudospins: a 
momentum-independent Rabi coupling (fla x ) and SMC 
(acr x sm(q x d), where q x is the quasimomentum and d 
the lattice period), with strengths of the same order. 
The coexistence of these two types of coupling leads to 
asymmetric FB band dispersions [38]. We investigate the 
FB band structures by measuring the quasimomentum 
of the BEC. The initial phase of the moving lattice 
plays a significant role in the Floquet dynamics [29], the 
effects of which are explored through a quantum quench 
induced dynamical coupling of the FB bands. Results 
are compared to theoretical predictions from a simple 
two-band model and from numerical simulations of the 
Gross-Pitaevskii (GP) equation. 

Results 

Experimental setup. To generate the s-p band SMC 
and FB band structures, we begin with a 87 Rb BEC 
composed of approximately 5 x 10 4 atoms confined in a 
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FIG. 1: Experimental setup and schematic lattice illus¬ 
tration. (a) Experimental arrangement. The crossed dipole 
trap beams propagate in the e x and e z direction. The static 
and moving lattice have overlapping beams propagating along 
e x +e y and —e x +e y . (b) Lattice potentials along the e x direc¬ 
tion. The lattice period d is identical for the static lattice Vo 
and the moving lattice V x . The initial offset between lattice 
sites of the static and moving lattice, Ax , is given by the ini¬ 
tial phase 0o between the two lattices. (c,d) Illustration of the 
multi-photon processes for the driven lattice system and the 
corresponding FB band structure in the first Brillouin zone. 
The static lattice induces a large energy gap (I) through a 
2-photon process and a small energy gap (II) through a 4- 
photon process. The moving lattice induces an energy gap 
when the s-band and p-band are coupled through (III). A 
smaller energy gap is produced by a combination of the static 
and moving lattice (IV). 


crossed dipole trap. A static lattice is generated by two 
perpendicular laser beams with wavelength A ~ 810 nm 
intersecting at the position of the BEC, as schematically 
shown in Fig. 1(a). The harmonic trap frequencies due to 
the envelope of the static lattice beams and the crossed 
dipole trap are (u: x ,ujy,uj z ) = 27r x (41,159,115) Hz, 
where e x points along the lattice, e y is the horizontal 
transverse direction, and e z is the vertical direction. A 
weak moving lattice with the same lattice period as the 
static lattice, d = tt/Icl where kL = v^tt/A, is then over¬ 
laid with the static lattice (Fig. 1(b)). The moving lat¬ 
tice beams are approximately 180 MHz detuned from the 
static lattice. A small frequency difference A w between 
the two moving lattice beams determines the velocity of 
the lattice according to vi at uce = A u /2kL- To induce s-p 
orbital band coupling, is chosen close to the energy 
gap E sp between the s- and p-bands of the static lattice 
at quasimomentum q x = 0. 

One outstanding feature of the coupling scheme em¬ 
ployed in these experiments is the asymmetry of the ef¬ 


fective s-p FB bands, which exhibit a local minimum lo¬ 
cated at a finite quasimomentum q x ^ 0. The direction 
in which the minimum is shifted away from q x = 0 is 
determined by the sign of (which determines the di¬ 
rection of motion of the moving lattice) and | A^ | — E sp 
(i.e. the detuning of the drive from the bandgap at 
q x = 0). Before describing experimental results and a 
formal derivation of the band structure using Floquet 
theory [29, 30], we lay the groundwork by presenting a 
multi-photon resonance picture that provides intuitive 
insights (Fig. l(c,d)). In this picture, one starts with 
the parabolic dispersion of a free atom in the absence of 
any external potentials. An optical lattice then induces 
2n-photon couplings (with n being an integer number) 
between points of the dispersion relation due to absorp¬ 
tion and stimulated emission processes. The couplings 
are centered around pairs of points that fulfill conserva¬ 
tion of energy and momentum. At these points, bandgaps 
open due to avoided crossings. Examples for possible cou¬ 
plings due to the static lattice (red arrows in Fig. 1(c)) 
and the moving lattice (blue arrows in Fig. 1(c)) and the 
associated bandgaps in the first Brillouin zone are shown 
in Fig. 1. Different coupling strengths lead to different 
sizes of bandgaps, which result in an asymmetric band 
structure. 

In another pictorial way, the Floquet band structure 
for the time-periodic system can be constructed by cre¬ 
ating multiple copies of the Bloch band structure of the 
static lattice that are offset in energy by | A w |. The mov¬ 
ing lattice couples the p-band and the shifted 5-band (la¬ 
belled by s' in Fig. 1(d)) at points where the shifted s- 
band intersects the unshifted p-band. The gaps opened 
by the coupling can formally be calculated using Floquet 
theory. 

Experimental measurements. Adiabatic loading of 
the BEC into an s-p FB band is achieved by first ramp¬ 
ing on the intensity of the static lattice, followed by adi- 
abatically ramping on the moving lattice intensity. In 
this procedure, the initial relative phase between the two 
lattices, 0o (Fig. 1(b)), becomes irrelevant and can effec¬ 
tively be set to zero. As we shall show in the context 
of Fig. 5, if the moving lattice is suddenly jumped on 
instead of adiabatically ramped on, this initial relative 
phase may manifest itself by drastically changing the dy¬ 
namics of the system [29]. 

Figure 2(a) shows the measured position, g min , of the 
band minimum for different driving frequencies, A w , af¬ 
ter adiabatically loading a BEC into a FB band. The 
driving frequencies are chosen such that hA^ lies in the 
gap at q x = 0 between the p-band (4.64 Er, where 
E R = h 2 k\/2m = h x 1749.5 Hz) and the d-band 
(5.44 Er). After adiabatically loading a BEC into a FB 
band, the lasers are switched off and the BEC is imaged 
after 14 ms time-of-flight (TOF). The positional shift of 
the BEC components is then used to determine the quasi¬ 
momentum. Each data point is an average over five iter- 











3 



A 0 [E r ] q x [k L ] 


FIG. 2: Effects of the driving frequency, (a) Band 
minimum qmin for the upper hybrid band vs. driving fre¬ 
quency Aw. The depth of the moving lattice is 1 Er. The 
filled circles are experimental measurements. The black line 
shows the theoretical prediction of a two-band model. The 
squares and stars are the results of numerical simulations of 
the Schrodinger equation and the GP equation, respectively, 
(b) Upper hybrid s-p FB band structure for different driving 
frequencies A w = 4.99 Er, 5.1 Er and 5.22 Er from top 
to bottom. The lowest (black) curve is the s orbital band 
without the presence of the driving field. 
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FIG. 3: Effects of the driving strength. Band minimum 
qmin vs. driving field strength V' for different driving frequen¬ 
cies of (a) |A W | = 2.92 Er and (c) |Ao,| = 5.21 Er. The red 
points are experimental data, the solid lines are the theoretical 
predictions from a two-band model. sgn(Auj) determines the 
direction of motion of the moving lattice. (b,d) Correspond¬ 
ing hybrid band structures for different driving field strengths 
Vx — 1.5 Er, 0.75 Er and 0 Er (outer to inner curves). 


ations of the measurement. A shift of the quasimomen¬ 
tum is detected that decreases with increasing driving 
frequency (Fig. 2(a)) as the coupling between the p-band 
and shifted s-band becomes weaker. The observed shift 
indicates a shift of the minimum of the upper hybrid band 
(Fig. 2(b)) into which the BEC is adiabatically loaded. 
The solid line in Fig. 2(a) shows qmin calculated from a 
simple two-band model (see below) and is in reasonable 
agreement with the data. The symbols are the results 
from real time simulation of the Schrodinger equation 
(squares) and the GP equation (stars) with finite nonlin¬ 
ear interaction strength [39] . We see that the interaction 
could modify the single-particle results. 

Figure 3 presents a complementary data set for which 
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FIG. 4: Heating of the Floquet system, (a) Number of 
atoms remaining after adiabatically loading a BEC into the 
FB band, normalized to initial atom number determined from 
independent experimental runs. The static lattice is ramped 
on to 5.47 Er in 200 ms. Then the moving lattice is ramped 
on to a depth of V' = 0.5 Er in 60 ms. The dips o, ft, and 
7 occur close to the Bloch bands p, d, and /. (b) Effective 
band structures for the lower hybrid band for data points 1, 
2, and 3 of panel (a), (c) TOF images taken at points 1, 2, 
and 3. 


the driving frequency is set to a constant value with 
| A^ | < E sp (Fig. 3a) or | A^ | > E sp (Fig. 3c) and the 
quasimomentum is determined for various depths of the 
moving lattice. The sign of A w determines the direction 
of motion of the moving lattice. For | A^ | < E sp the BEC 
resides in the lower hybrid s-p FB band (Fig. 3(b)) while 
for lA^I > E sp it is in the upper hybrid band (Fig. 3(d)). 
This leads to a shift of the quasimomentum into oppo¬ 
site directions for the two cases. For a given driving 
frequency, the coupling of the two bands is stronger for 
larger driving field strength (i.e. larger depth of the mov¬ 
ing lattice) so that the BEC is shifted to a larger absolute 
value of quasimomentum. 

Floquet systems such as the one in our experiment 
are described by quasienergy bands. They do not have 
a thermodynamic ground state, and in the presence of 
many-body interactions their stability can be affected by 
a variety of factors [40-42]. Experimentally, we study 
the stability of the system by determining the number of 
condensed atoms left after the static and the moving lat¬ 
tices are successively and adiabatically ramped on. TOF 
imaging reveals atom loss and heating of the BEC as 
shown in Fig. 4. The dips a , /?, and 7 in Fig. 4(a) 
occur when the driving frequency is chosen such that it 
leads to a coupling close to the Bloch bands p, d and / 
of the static lattice at q x = 0 respectively. The lower 
hybrid band structure, in which the BEC mainly resides, 
for points 1, 2, and 3 and the corresponding TOF im¬ 
ages are shown in panels (b) and (c). Resonance induced 
collective excitations and modulational instabilities can 
play a role for the observed losses [43] . 

Minimal two-band model. The dynamics of the BEC 
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are governed by the full time-dependent GP equation, 
ink'll)(r,t) = [H 0 (t) + Vtrap + V in t]-ip{r,t ) where V trap and 
Vi n t are the external trapping potential and the mean- 
field interaction, respectively. Ho(t) is the single-particle 
Hamiltonian, 

Hq (t) = ^—+ Vo cos 2 (k L x) + Vcos 2 {k L x + <j> 0 - 

(!) 

where the second and the third terms describe the static 
and moving optical lattices, respectively, and 0o is the 
initial relative phase between the two sets of lattices. 

When the static lattice depth Vo is large and when 
|A W | is close to the energy gap E sp , higher orbital bands 
are not significantly populated in the driven process and 
the system is well described by a simple two-band tight- 
binding model [38]. Following the standard procedure 
in Floquet theory, we obtain the effective single-particle 
Hamiltonian 


where 


rreff 

^0 


€s(Qx) A sp \ 

A *sp £p(Qx) - |AJ ) ’ 


( 2 ) 


A sp = —i[Q - a sin (q x d) + ft cos(q x d)]e (3) 

is the coupling between s- and p-orbital bands that is 
induced by the moving lattice potential for > 0 [44], 
and e s and e p are the energy dispersions for the uncou¬ 
pled orbital bands. The three coupling coefficients Q, 
a and (3 are given by Q = ^-(si\sm(2kLx)\pi), a = 
V ( Si | cos(2fc L a;)|pi + i) and /3 = ^-(si\sm(2k L x)\p i+1 ), 
where |s^) and \pi) are the maximally localized Wannier 
orbital states in the i-th site. Q is the coupling between 
s- and p- orbital states in the same lattice site, while a 
and (3 are the couplings between s- and p-orbital states 
of nearest neighbouring sites. SMC between s-p band 
psuedospins is represented by <asin (q x d)a x . 

This derivation shows that the inversion symmetry of 
FB band structure is broken due to the coexistence of 
couplings of different parities. When the moving lat¬ 
tice depth is adiabatically ramped on, the quasimomen¬ 
tum of the BEC gradually shifts away from q x = 0 in 
a definite direction following the hybrid band minimum. 
This is quite different from previous shaken lattice ex¬ 
periments [17] where the inversion symmetry of the band 
was preserved and the BEC could spontaneously choose 
either side of q x = 0 as its ground state. In that case, 
the BEC needed to be accelerated to break the inver¬ 
sion symmetry. In our scheme, the position of the true 
minimum is uniquely determined by the moving velocity 
direction, moving lattice depth, and driving frequency. 

This minimal two-band model captures the essential 
physics of the driven lattices as we have seen through 
the comparison of experimental measurements and theo¬ 
retical values (see Figs. 2 and 3), demonstrating the ob¬ 
servation of SMC between s-p band pseudospins. How¬ 
ever, this model may deviate from the experiment when 
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FIG. 5: Quench dynamics after suddenly jumping on 
the coupling between the s and p band, (a, b) Nor¬ 
malized occupation of the momentum component —2 hk (blue 
points with solid error bars) and +2 hk (red points with dashed 
error bars) in (a) and 0 hk in (b). The dots are the average of 
ten experimental measurements for each time. The error bars 
indicate the spread of the experimental data. The shaded 
areas are the results of numerical GP simulations calculated 
for a homogeneous distribution of different initial phases <3 o- 
The black curve represents the calculation for phase (f>o = 0 

(c) Bandstructure plot. Jumping on the moving lattice places 
the BEC (black ellipse) into the gap between two FB bands. 

(d) Experimental images taken 0.5 ms after the quench for the 
top two images and 0.8 ms after the quench for the bottom 
image. Dashed squares indicate the areas used for counting 
the atom number in the —2 hk (left square), 0 (middle square) 
and +2 hk component (right square). The sum of the atoms 
in all three boxes is used for the normalization of the experi¬ 
mental data in panels (a) and (b). 


the modulated dynamics involve additional orbital bands 
or when the nonlinear interaction is strong such that the 
single-particle band structure will be renormalized by the 
interaction term. 

Quench dynamics. Since a Floquet system is gener¬ 
ated by a time-periodic Hamiltonian, an important ques¬ 
tion concerns the role of the initial phase of the driv¬ 
ing field [29]. For the system considered in this work, 
this phase determines the relative positions between the 
moving and static lattice sites. Though the relative phase 
does not change the effective band structure (Eq. 2), and 
thus the time-averaged dynamics, it can play a crucial 
role in the micromotion of the BEC. To demonstrate the 
effect of the initial relative phase, we study the oscilla¬ 
tions in the population of the momentum components 
k x = 0,^2/^ after a quantum quench. Figure 5 (a-b) 
present such quench dynamics after adiabatically ramp¬ 
ing on the static lattice to 5.47 Er followed by a sudden 
jump on of the moving lattice to V x = 1 Er with an 
on-resonant driving frequency || = E sp (Fig. 5 (c)). 
We focus on the evolution during the first 3 ms, during 
which the BEC mainly stays at q x = 0 without signif- 
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icant dipole motion in the hybrid bands. The symbols 
in Fig. 5 (a-b) are experimental data averaged over ten 
measurements for each time step. There is significant 
spread in the data for each time step, as indicated by 
the vertical error bars. This spread is due to the ini¬ 
tial phase 0o between the static and the moving lattice, 
which is uncontrolled in the experiment, such that each 
iteration realizes a case with a different, random 0o- The 
shaded areas represent the result of numerical GP simu¬ 
lations for a homogeneous spread of relative phases. The 
experimental error bars are in reasonable agreement with 
the expectation based on these numerics. The numerics 
reveal that for a fixed initial phase there are two oscilla¬ 
tion periods of different timescales (Fig. 5(b)). The fast 
oscillation (of period T ~ 0.1ms) corresponds to the mi¬ 
cromotion of particles under the high-frequency periodic 
driving, whereas the slow oscillation (T ~ 1.75ms) cor¬ 
responds to the time-averaged effective Rabi oscillations 
between the two hybrid FB bands. For longer holding 
time, the periodicity is slightly broken due to a small 
dipole motion. 

Discussion. We have realized and characterized a new 
kind of SMC with lattice bands as pseudospins. This not 
only provides a powerful tool to control orbital states 
with a driving field, but also enriches the study of novel 
quantum matter using hybrid orbital bands. There are 
many directions that can be taken along this route, e.g., 
the engineering of similar SMC in higher dimensional sys¬ 
tems involving different orbital bands, and quantitative 
analysis and measurements of the effects of strong inter¬ 
actions on the effective bands. The realization of similar 
SMC for fermionic atoms such as 6 Li and 40 K with tun¬ 
able interactions may open the door for exploring exotic 
quantum matters such as FFLO superfluids and Majo- 
rana fermions. 
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